LHCb Masterclass : Measuring the D lifetime at the LHC 



1 Introduction 

The Large Hadron Collider (LHC) is not only a tool with which to search for exotic new particles, 
but also a factory of particles whose existence is in no doubt but whose precise properties are not 
yet known well enough. One example are charmed particles, i.e. particles which contain a charm 
(also known as c) quark, which were first discovered more than 30 years ago. Approximately one in 
every ten LHC interactions produces a charmed particle, and the LHCb experiment at the LHC has 
already recorded over one billion events containing charmed particles. 

Even though the LHC produces a lot of charmed particles, it produces a much greater number of 
other particles which can be mistaken for the signal ( "backgrounds" ). In order to extract information 
from such large signal samples, it is necessary to achieve excellent control over these backgrounds. 
Today you will be performing exercises using real collisions recorded by the LHCb experiment during 
2011 datataking, which contain both signal and background particles. This set of exercises is 
designed to teach you how to 

1. Use an event display of the proton-proton collisions inside the LHCb detector to search for 
charmed particles and separate this signal from backgrounds. 

2. Fit functional forms for the signal and background to the data in order to measure the number 
of signal events in the data sample and their purity (defined as the fraction of signal events in 
the total sample). 

3. Obtain the distribution of signal events in a given variable by taking the combined distribution of 
events in the data sample (which contains both signal and background events) and subtracting 
the background distribution. The result of the fit in the previous step is used to find a sample 
of pure background events for subtraction, and to compute from the signal yield and purity the 
appropriate amount of background which should be subtracted. 

4. The signal you will be looking at decays exponentially with time, analogously to a radioactive 
isotope. You can now use the sample of events passing the previous step to measure the 
"lifetime" of the signal particle. The lifetime is defined as the time taken for (e — l)/e of the 
signal events to decay, where e « 2.718 is the base of the natural logarithm. It is analogous 
to the concept of half-life in radioactive decay. 

The data sample used for this exercise consists of candidates for a type of charmed particle known 
as a D° particle found in a sample of randomly collected LHC interactions during 2011 datataking. 
A D° particle consists of a charm quark and an up antiquark. The particles are measured decaying 
in the mode D° — > K~ir + , where the final state particles are a kaon (K~) consisting of a strange 
quark and an up antiquark, and a pion which consists of a down antiquark and an up quark. 
These particles have lifetime which are long enough that, for the purpose of this exercise, they are 
stable within the LHCb detector. The particles have been preselected using loose criteria so that 
you begin the exercise with a visible signal, but background events are also present. 

Before discussing the exercise further, it is worth taking the time to familiarize yourselves with 
the LHCb detector, shown in Fig. 1. It is a forward spectrometer covering the angular range between 
0.7° and 15° relative to the LHC beamline. The beamline is located at y = 0 in the picture and runs 
in the z direction. In what follows "transverse" means transverse to the LHC beamline, "upstream" 
is towards smaller values of z and "downstream" is towards greater values of z. The detector 
includes a high-precision system for tracking charged particles consisting of a silicon-strip detector 
surrounding the proton-proton interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip detectors 
and straw drift tubes placed downstream. Charged particles leave straight line tracks in the detector 
surrounding the interaction region, where there is no magnetic field, and are subsequently bent by 
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Figure 1: The LHCb detector. The Z axis is the direction of the LHC bcamlinc. 




Figure 2: The D° particle secondary vertex ("D") displaced from the pp interaction (yellow star-like blob). 

the magnet before leaving tracks in the downstream tracking station. Their momentum and charge 
can be deduced from the curvature induced by this magnetic field. The combined tracking system 
has a momentum resolution Ap/p that varies from 0.4% at 5 GeV/c to 0.6% at 100 GeV/c, an 
impact parameter 1 resolution of 20 /zm for tracks with high transverse momentum, and a lifetime 
resolution of 50 fs. 

A distinctive feature of D° particles is their measurably long lifetime (the object of this exercise). 
From a practical point of view, this means that they form a "secondary" vertex which is displaced 
from the "primary" proton-proton interaction. This is illustrated in Fig. 2. This fact, together with 
their abundant production rate mentioned earlier, allows the D° signals to be well separated from 
the background of the underlying event, most of which consists of random combinations of particles 
produced in the proton-proton collision. 



1 Impact parameter is the transverse distance of closest approach between a particle trajectory and a vertex, most commonly 
the primary proton-proton interaction vertex. 
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Figure 3: Three views of one event, from top to bottom : y-z, y-x, and x-z. 
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The data is provided together with a graphical user interface (GUI) which you can use to access 
it - go ahead and execute it. You can access online instructions from the GUI using the "help" 
button. The exercise consists of two parts : the first is an event display which you can use to search 
for the displaced vertices of D° particles inside the LHCb detector, the second is a fitter which you 
can use to separate signal from background and measure the lifetime of the D° particles. 

2 Event display exercise 

The aim of the event display exercise is to locate displaced vertices belonging to D° particles in the 
vertex detector of the LHCb experiment. When you launch the exercise and load an event, you will 
see an image of the LHCb detector and particle trajectories ("tracks") inside it. These tracks are 
colour coded, and a legend at the bottom of the GUI tells you which colour corresponds to which 
kind of particle. 

In order to make identifying vertices easier, you can view an event in three different two- 
dimensional projections : y-z, y-x, and x-z, show for one event in Fig. 3. Different events will 
be clearer in different projections, so feel free to experiment with all three! Displaced vertices ap- 
pear as a pair of intersecting tracks, far away from the other tracks in the event. When you click 
on a particle, you will see its information, including mass and momentum, in the Particle Info box. 

A D° particle decays into a kaon and a pion, so you will need to find a displaced vertex where a 
kaon track intersects with a pion track. Once you find a track which you think is part of the displaced 
vertex, you can save it using the Save Particle button. Once you have saved two particles, you can 
compute their mass by clicking on the Calculate button. If you think this combination has a mass 
compatible with that of the D° particle, click on Add to save it : by saving a combination for each 
event, you will build up a histogram of the masses of the displaced vertices in the different events. 

Remember that you are looking at real data so it contains both signal and background, and the 
detector has a finite resolution, so not all displaced vertices will have exactly the D° mass (even the 
signal ones). They should, however, be within the range 1816-1914 MeV (this range is around 3% 
each way around the true D° mass). If you try to save a combination which is too far away from 
the real D° mass the exercise will warn you that you have not found the correct displaced vertex 
pair and won't let you save it. 

If you are not able to find the displaced vertex for an event after a few minutes, move on to the 
next event and come back to the one which was giving you trouble if you have time at the end of 
the exercise. Once you have looked at all events, you can examine your mass histogram by clicking 
the Draw button. Discuss the shape with a moderator. 

At the end, please remember to click on Save Histogram in order to save your histogram to 
disk! It will then be combined by the moderators with the histograms of the other students and you 
should discuss your results as a group. 

3 Fitting exercise 

Before describing the fitting part of the exercise, it will be useful to list the variables involved in this 
exercise : 

• DO mass : this is the invariant mass of the D° particle. The signal can be seen as a peaking 
structure rising above a flat background. The range of masses relevant for this analysis is 
1816-1914 MeV. The signal shape is described by the Gaussian (also known as "normal") 
distribution. The centre ("mean") of this distribution is the mass of the D° particle, while the 
width represents the experimental resolution of the detector. 

• DO TAU : this is the distribution of decay times of the D° candidates. The signal is described 
by a single exponential whose slope is the D° lifetime (the object of the last exercise), while 
the background concentrates at short decay times. 

• DO IP : this is the D° distance of closest approach ("impact parameter") with respect to the 
proton-proton interaction in the event in millimetres. The smaller the impact parameter, the 
more likely it is that the D° actually came from that primary interaction. In order to simplify 
the drawing, we actually plot and cut on the logarithm (base 10) of this quantity in the exercise. 
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• DO PT : this is the momentum of the D° transverse to the LHC beamline. 
Now let's proceed through the exercise itself. 

3.1 Exercise 1 : fitting the mass distribution and obtaining signal variable distri- 
butions 

The object of this exercise is to fit the distribution of the DO mass variable, and extract the signal 
yield and purity. 

1. Click on the Plot DO mass button to plot the overall mass distribution. You will see a peak 
(signal) on top of a flat distribution (background). The peak should be described by a Gaussian 
function, whose mean corresponds to the mass of the D° and whose width (a) is determined 
by the experimental resolution of the LHCb detector. 

2. Click on Fit mass distribution to fit this distribution using a Gaussian function for the signal 
and a linear function for the background. 

3. Look at the fitted mass distribution. You can split it into three regions : the signal region and 
two background-only "sidebands" : one above the signal (the upper sideband) and one below 
the signal (the lower sideband). A Gaussian distribution contains 99.7% of its events within 
three standard deviations of the mean, so this "three sigma" region around the mean is usually 
the definition of the signal region. 

4. Use the slider labelled Sig range to define the beginning and end of the signal region. All 
events not falling into the signal region will be said to fall into the background region. 

5. You can now use the definitions of the signal and background regions in the mass variable 
to determine the signal and background distributions in other variables. Click on the button 
labelled Plot distributions. You will see the signal (blue) and background (red) distributions 
for the other three variables plotted next to the mass distribution. You should discuss the 
exercise with an instructor at this point. 

3.2 Exercise 2 : measuring the D° lifetime 

The object of this exercise is to use the signal sample which you obtained in the previous step to 
measure the lifetime of the D° particle. This is the same quantity as the half-life of a radioactive 
particle : the D° decays according to an exponential distribution, and if this exponential is fitted to 
a distribution of the D° decay times, the slope of this exponential is the reciprocal of the lifetime 
of the D°. 

1. Fit the lifetime of the D°. 

2. Compare the slope of this exponential to the reciprocal of the D° lifetime given by the Particle 
Data Group. Talk to an instructor about how well these agree with each other. 

3. In addition to statistical uncertainties, measurements can suffer from systematic uncertainties, 
caused by a miscalibrated apparatus or an incorrect modelling of the backgrounds. One basic 
technique for estimating these is to repeat the measurement while changing the criteria used 
to select signal events. If the result changes significantly when changing the criteria, we know 
that there is something wrong! 

4. Repeat your fit for the lifetime of the D° while varying the maximum allowed logarithm of the 
D° impact parameter. The allowed values range from —4.0 to 1.5 in the original fit. Move this 
upper value from 1.5 to —1.9 in steps of 0.20, and refit the D° lifetime at each point, saving 
the results as you go along. 

5. Plot the histogram showing the fitted value of the D° lifetime as a function of the upper cut 
on the logarithm of the impact parameter. Discuss the shape, and what it tells you about the 
D° lifetime, with an instructor. 

6. What other sources of systematic uncertainty might we need to consider when making a lifetime 
measurement? 
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